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We show that the gauge-invariant transverse-momentum-dependent (TMD) quark distribution 
function can be expressed as a sum of all higher-twist collinear parton matrix elements in terms 
of a transport operator. From such a general expression, we derive the nuclear broadening of the 
transverse momentum distribution. Under the maximal two-gluon correlation approximation, in 
which all higher-twist nuclear multiple-parton correlations with the leading nuclear enhancement 
are given by products of twist-two nucleon parton distributions, we find the nuclear transverse 
momentum distribution as a convolution of a Gaussian distribution and the nucleon TMD quark 
distribution. The width of the Gaussian, or the mean total transverse momentum broadening 
squared, is given by the path integral of the quark transport parameter qF which can also be 
expressed in a gauge invariant form and is given by the gluon distribution density in the nuclear 
medium. We further show that contributions from higher-twist nucleon gluon distributions can 
be resummed under the extended adjoint two-gluon correlation approximation and the nuclear 
transverse momentum distribution can be expressed in terms of a transverse scale dependent quark 
transport parameter or gluon distribution density. We extend the study to hot medium and compare 
to dipole model approximation and N = 4 Supersymmetric Yang-Mills (SYM) theory in the strong 
coupling limit. We find that multiple gluon correlations become important in the strongly coupled 
system such as M = 4 SYM plasma. 

PACS numbers: 13.85.Hd, 25.75.Bh,11.80.La 



I. INTRODUCTION 



The success of perturbative QCD (pQCD) in describing hard processes in hadronic interactions relies on the fac- 
torization theorem [l[ that separates the coherent long distance interaction between projectile and target from the 
incoherent short distance interactions. The physical observables such as cross sections of deeply inelastic scattering 
(DIS) and Drell-Yan (DY) dilepton production can be expressed as a convolution of hard partonic scattering cross 
sections, parton distribution functions and parton fragmentation functions. The hard partonic parts are calculable in 
a perturbative expansion in the strong coupling constant a s (Q 2 ) which becomes small large momentum scale Q 2 of 
the hard processes 0, 9- Though the parton distribution and fragmentation functions are not calculable in pQCD 
since they involve long distance interaction, they are universal and independent of the specific partonic hard processes. 
Therefore, they can be measured in one hard process and then applied to another, therein lies the predictable power 
of pQCD. 

The most practiced factorization scheme is collinear factorization in which one integrates out the transverse mo- 
mentum of the initial (final) parton up to a factorization scale and the final observables will only depend on the 
transverse-momentum-integrated or collinear factorized parton distribution (fragmentation) functions. Such a proof 
of factorization has also been extended to semi-inclusive processes [U, [1J that involve finite transverse momentum of 
the final hadron or dilepton with the introduction of transverse-momentum-dependent (TMD) parton distribution 
and fragmentation functions. The final observables can be expressed as a convolution of collinear hard parts (setting 
the initial parton transverse momenta to zero) and TMD parton distribution and fragmentation functions |6fl. Such 
TMD parton distribution and fragmentation functions are important for the study of hadronic interactions with singly 
or doubly polarized beams, such as single-spin asymmetry in semi-inclusive processes in DIS (SIDIS) @, 1, 1 and 
proton-proton scattering. 

In the proof of factorization [lOj], in DIS off a nucleon or nucleus target for example, one important step is to 
cikonalize all soft interactions between the struck quark and the target remnant as shown in Fig. [TJ The summation 
of these soft gluon interactions gives rise to a definition of TMD quark distribution function in a nucleus, 

/,W±) = / ^-^e^y'-^^iA | ^(0,0j^£ T MD(0,y)^-,^) | A), (1) 
where, 

£tmd(0, y) = d (-oo, 0; ± )£^(-oo; y x , O x )£ii (-oo, y-;yx) (2) 



XnP 



Ap 



Vn 




%2P Xip 



xp 



2/2 yi 



FIG. 1: Multiple soft gluon interaction between the struck quark and the remnant of the target nucleus in DIS. 
is the complete gauge link in TMD quark distribution function that contains both the transverse 
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d£± ■ A_l(-oo,£_l) 



(3) 



£\\(-oo,y ;y±_) = Pexp 



-'9 I d£ ,V±) 
>y~ 



(4) 



The above gauge links for quark propagation are defined in the fundamental representation, = A^T a (a = 1 — 
8). Note that in our convention, the quark is produced at (y~,y±) [ or (0,0j_)] and propagates toward (— oo,y±) 
[(— oo,0j_)]. The path ordering along the light-cone is then defined from y~ to — oo. 

Even though the parton distribution in Eq. ([T]) describes the probability to find a quark with momentum fraction x 
and transverse momentum k±_ in a nucleus, it also contains information about final state interaction with the target 
remnant encoded through the gauge links. These gauge links are not only crucial to ensure the gauge invariance of 
the TMD parton distribution functions in both light-cone and covariant gauge but also lead to physical consequences 
such as single-spin asymmetry in semi- inclusive DIS and Drell-Yan process [8. IToL [ill] . For DIS off a nucleus target, 
they should also contain information about transverse momentum broadening of the struck quark due to multiple 
scattering inside the nucleus. The main purpose of this work is to study nuclear transverse momentum broadening 
in DIS from the TMD quark distribution functions and extend the result to the case of quark or jet propagation in a 
thermal medium as in high-energy heavy-ion collisions. 

In the study of nuclear matter and quark-gluon plasma in high energy lepton- nucleus, hadron- nucleus and nucleus- 
nucleus collisions, jet transverse momentum broadening plays a crucial role in unravelling the medium properties 
through modification of the final jet or hadron spect ra ( jet quenching) due to final state interaction between the 
energetic partons and the nuclear or hot medium fl2l . fl3| . Current phenomenological studies of experimental data 
EE Ell E3, El] on jet quenching rely on pQCD calculations of the parton energy loss or modification of the parton 
fragmentation functions due to gluon radiation induced by multiple scattering during the parton propagation in the 
medium. One important parameter that controls parton energy loss or medium modification of the jet fragmentation 
function is the jet transport parameter q or transverse momentum broadening squared per unit of propagation length 
(l9l I20I I2H I22I I23I [U [2||. Therefore, calculation and measurement of the jet transport parameter is an important 
step toward understanding the intrinsic properties of the QCD medium. 

There are many calculations of the jet transverse momentum broadening. Early studies dealt with the mean average 
transport parameter [2ll . |26| under one particular gauge without apparent guaranty of gauge invariance. Transverse 
momentum broadening in Drell-Yan pair production in proton-nucleus scattering (27j and a recent calculation of 
nuclear transverse momentum broadening distribution in DIS J28J were obtained by a direct summation of multiple 
scattering in covariant gauge, again without apparent gauge invariance in the final result. One approach to the 
parton propagation in medium with a gauge (mostly) invariant framework is the Wilson line formulation of multiple 
parton scattering [23| that resembles the longitudinal gauge link in Eq. ([1]), but again in the covariant gauge. Since 
resummation of all soft gluon interactions is crucial for the gauge invariant form of the TMD parton distribution 
function in Eq.fl]), it must contain the nuclear transverse momentum broadening due to multiple parton scattering 
in nuclei. One should be able to derive a completely gauge invariant form of the nuclear transverse momentum 
broadening from the nuclear TMD quark distribution function. This is what we will prove in this paper. We will 
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derive the nuclear broadening of the transverse momentum distribution simply from the gauge invariant form of 
TMD quark distribution function in Eq.|T]) with a given nuclear distribution function. The broadened distribution 
will have a Gaussian form as found in earlier studies, considering only two-gluon field correlations in the nucleons. 
However, the broadening parameter in our derivation will have an explicit gauge invariant form. We will also show 
that contributions from higher-twist multi-gluon correlations in a nucleon can be rcsummed to give a nuclear TMD 
quark distribution that depends on the transverse scale dependent gluon distribution density inside the nucleus. 

We will briefly summarize our main results here. To calculate the nuclear transverse momentum broadening 
including all higher-twist nuclear parton matrix elements, we first express the gauge-invariant TMD quark distribution 
function in a nucleus as a sum of collinear higher-twist nuclear parton matrix elements which can be exponentiated 
to give, 

ff(x B ,k ± ) = J ^-e^ + y~(A | ^(0,0 ± )^£ | |(0,y-;0 ± ) e # -^^ 3 ->^^(^,0 ± ) | A)6^(k ± ), (5) 
where the transport operator W±(y~ ,y±) is defined as 



W ± (y-,y ± ) = iD ± (y-,y* x )+g d£~ C\{£,y- ,y±_)F +± (C ,y±.)C { \(C ,y~;y±), (6) 

J — oo 

and D±(y~ ,y±) = d± + igA±(y~ ,y±) is the covariant derivative. The transport operator is supposed to operate 
on both the quark field and the gluon fields within itself and the transverse coordinate is set to y*±_ — 0± after the 
operation. 

With a maximal two-gluon correlation approximation, the high-twist multi-parton correlations in a large and weakly 
bound nucleus can be expressed as products, which have the maximum nuclear size (or medium length) dependence, 
of twist-two nucleon gluon distribution functions. The leading contribution to the nuclear TMD quark distribution is 
shown to have a simple form, 



fq(x B , kx) = Aexp 



d£ N q F {£, N ) 



H 7tA 2 f 



with A2F = J d£,N^ F (^ N ) as the total transverse momentum broadening squared. Such a distribution also satisfies 
a 2-D diffusion equation in transverse momentum k± with the diffusion constant given by qf, the so-called (twist-2) 
quark transport parameter, 

MU) = ^Jp4(6v) J ^(^(^iTr^^O^Kii^rsOi^+^r-O^Kj^r;^)] I iV), (8) 

which is related to the twist-two collinear nucleon gluon distribution function. It is also the mean transverse momentum 
broadening squared per unit length and is proportional to gluon distribution density inside the nucleus. Here P^{£,n) 
is the single nucleon density inside the nucleus and p + is the longitudinal momentum per nucleon. Inclusion of 
higher-twist nucleon gluon matrix elements under an extended adjoint two-gluon correlation approximaton will give 
rise to a similar nuclear transverse momentum distribution with an effective transverse-scale-dependent = — ) 
transport parameter, 

fr(6r,£) = ^/4(6v) J H+(- 2 )^ I Tr [F + A0,0 ± )C TM u(0,OF + ^r,^K MD (0,O] I N), (9) 

which is related to the transverse-scale-dependent gluon distribution function in a nucleon. Such transverse-scale- 
dependent transport parameter constituents power corrections to the Gaussian form of the final nuclear broadening 
distribution. 

The rest of this paper is organized as follows. In the next section, we will first formulate the nuclear TMD quark 
distribution in terms of collinear nuclear high-twist parton matrix elements in the light-cone gauge and then derive 
the nuclear transverse momentum broadening distribution under the maximal two-gluon correlation approximation. 
We then include the higher-twist nucleon quark matrix elements to consider the effect of nucleon intrinsic transverse 
momentum. The derivation of nuclear broadening of the TMD quark distribution is generalized to any arbitrary 
gauge with the final results expressed in an explicit gauge invariant form. Effects of higher-twist gluon distribution 
functions are also discussed and are shown to lead to a transverse-scale-dependent quark transport parameter. In 
Sec. IIII1 we will extend the results on nuclear transverse momentum broadening to the case of quark propagation in 
hot medium and compare our results under the maximal two-gluon approximation to the dipole model approximation 
in the Wilson line approach to multiple parton scattering. We will also compare to the result from J\f = 4 super- 
symmetric Yang-Mills (SYM) calculation 30] of the transverse momentum broadening and discuss the importance of 
multiple gluon correlation in a strongly coupled system. Finally we will give a summary in Sec. HVi 
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II. NUCLEAR TRANSVERSE MOMENTUM BROADENING 



We will consider DIS off a large nucleus as depicted in Fig. [TJ In the infinite momentum frame, the nucleus has 
a longitudinal momentum p = [p + ,0, 0_J per nucleon and a quark with fractional longitudinal momentum ig = 
Q 2 /2p + q~ is knocked out by a virtual photon with four-momentum q — [— Q 2 /2q~, q~ , OjJ. The differential cross 
section for e~ {h) + A(Ap) — > e~(l 2 ) + q(k) + X can be written as, 

w W da - "em 1 r^n j \j? rif y. 



where s = (p + l\) 2 and alpha rm EM is the fine structure constant in electrodynamics. The SIDIS hadronic tensor 
Wpv is defined as, 

Ek^ = Y,( A \Mm,X)(k,X\M0)\A)2TrS i (p + q-k-p x ), (11) 
x 

and Jfj.(y) = l >p(y)'YiStp(y) IS t ne hadronic electromagnetic current. The leptonic tensor L^ v is defined as usual and is 
given by, 

L^(l u l 2 ) = ~Tr[Z/ 7/ ^ 7/1 ] = 2[ZjV + 15V " (h • h)gT\- (12) 

The struck quark carrying large negative longitudinal momentum q~ will suffer multiple soft scattering with the 
rest of the nucleus before hadronization into hadrons. We assume that the virtuality of the photon Q 2 is very large 
and consider now only the lowest order of the hard partonic part. The soft interaction as shown in Fig. [T] can be 
resummed and the final leading twist SIDIS tensor 0] 

^ = HjS(x B p,q)ff(x B ,k ± ), (13) 

can be factorized as the product of the lowest hard partonic part H^J (xbP, <?) , 

H$(x B p, q) = ^Trthr(x B f + i)lA^~ q (14) 

and the nuclear TMD quark distribution function fg(x B ,k±) as defined in Eq. |lj. Since the final state interaction 
are already included in the nuclear TMD quark distribution function as the gauge links, one should be able to derive 
the nuclear transverse momentum broadening. 



A. Nuclear TMD quark distribution function in light-cone gauge 

One important feature of the complete and gauge invariant nuclear TMD quark distribution function in Eq. ([T]) is 
the transverse gauge link £j_(— oo; y±, 0±) which depends on the transverse gauge potential A±(— oo, at — oo along 
the light-cone. Without it, the gauge links in the TMD quark distribution would completely vanish in the light-cone 
gauge and one will be misled to assume that the final state interactions become absent. Furthermore the TMD quark 
distribution without the transverse gauge link is no longer gauge invariant under residual gauge transformation since 
the transverse gauge potential at infinity does not vanish in the light-cone gauge [3l| and is closely related to the 
singularity of the gluon propagator in the light-cone gauge which has to be properly regularized. Therefore, the 
effects of final state interaction are actually encoded in the transverse gauge link in the light-cone gauge and cannot 
be casually discarded. For this reason, we will first derive the nuclear transverse mometnum broadening in light-cone 
gauge and repeat the derivation later in an arbitrary gauge. 

Let us first consider the nuclear TMD quark distribution function in the light-cone gauge A+ = 0. In this gauge all 
the longitudinal gauge links vanish and we are left with only the transverse gauge link in the TMD quark distribution 
function. We first insert a S function into the TMD quark distribution function, 



f A ( X ,k ± ) = J d 2 t ± tf(x,i ± )6W(k ± -i ± ). 



(15) 
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Using a Taylor expansion of the <5-function, 

the quark transverse momentum distribution can be written as 
dy~ d 2 y± 



(16) 



2ir (2tt) : 



:d 2 £±e lxp v ~ tl± - 



x(A | ^(0,0 ± )^4(-oo;^,0 ± My-,^) I A) e - e ^^S^(k ± ) 
^^d 2 t ± e^ + y-- ir ^(A | t>(0,d ± )l-e> B y^^c{(-™; yx, Ox)V>0T, £i) I ^ (2) (fc±), (17) 

after partial integration in the transverse coordinate yx- Since both the quark field and the transverse gauge link 
depend on y±, we have 



id yx C' x (-oo;y±,fix) = £^(-00; y ± , 0±)[-gA±(-oa, y±_) + id y± ] 



C]_(-oo;yj_,0±) 



iD±(y ,y±)+9 d£ ^+x(£ ,y ± ) 



(18) 



where Dj_(y ,y±) = d y± +igA±(y ,y±) is the covariant derivative. In the last step of the above equation we used 
the following identity, 



Ai(-oo-,y x ) = A±(y-,y x )- / dCd+A ± (C ,y±), (19) 

J — 00 

and d+A± — i*+x in the light-cone gauge. Completing the integration over the transverse momentum i\_ in Eq. \YI\ 
will now produce a S- function 5^ 2 ^(yx) which will set transverse coordinate y± = Ox at which the transverse gauge 
link will disappear. We have then, 



f q A (x,k ± ) = J d JL-e^v- {A 1 ^ ( o)XL e ^xGr)-v^ (jr) 1 A )S^(k±). 

Here we define the transport operator W±(y~,y±) in the light-cone gauge as 



(20) 



W ± {y-,y ± ) = iD ± (y-,y ± )+g / dTF+±(C,&±)- 

J —OO 

Note that in the light-cone gauge, the transport operator is translational invariant, 

ti r ±(y-,y±) = Wx(o,yx), 



(21) 



(22) 



along the light-cone. For brevity in notation we will suppress the transverse coordinates whenever they are set to zero 
in the field operators, 



O(y-,0 x )=O(y-). 



(23) 



Eq. (|20p is a general result for the transverse momentum quark distribution function in the light-cone gauge. In the 
case of a large nucleus target, we can make further simplifications under the assumption of a weakly bound nucleus. 
We first expand the exponential factor in Eq. (|20p in power of the transport operator W±(0). The expectation value 
of any odd power of the operator under any unpolarized nuclear state should vanish under the parity invariance. We 
therefore are left only with the even-power terms of the expansion, 



1 - ^ + r -> — 1 2 ™ 

M 2n = j2ny_i A \ ^(Q)-y [w±(y-) • V fci J i>{y~) I A). 



(24) 



We will neglect the covariant derivative first in the transport operator and consider first the twist-four nuclear 
matrix elements (n — 1), 



d£(A I m\F +i (^)F +j (^(y-) I A). 



(25) 
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Because a nucleus consists of nucleons which are color singlet states, the quark and gluon fields could either be all 
attached to a single nucleon or to two separate nucleons. In the first case, all four parton fields in the above correlation 
matrix elements are confined to the size of a nucleon y~ , , ~ r N- On the other hand, if quark and gluon fields are 
confined to two separate nucleons, y~ , = |£f — £<f| ~ rjv, the overall position of the gluon field £Z = (£f + £^~)/2 
will follow the second nucleon and are only confined to the size of the nucleus Ra- Therefore, the quark-gluon 
correlation function in this case will have a nuclear enhancement of the order Ra/tn ~ A 1 / 3 as compared to the first 
case where both quark and gluon fields are confined to a single nucleon. As a two-parton correlation approximation 
for a large nucleus target we will only keep the matrix elements with the nuclear enhancement. We will also neglect 
the correlation between different nucleons and assume the large nucleus as a weakly bound. The leading contribution 
to the above quark-gluon correlation function will be then, 

M 2 * A J dy-e^ + y~(N | ^(O)^(jT) | N) 9 ^ J d£f £ d£ ^-«Tr [F +i (^)F +j (^)])) A . (26) 

If we further assume the large and weakly bound nucleus as a homogenous system of nucleons, 

((Tr [F +i (^)F +j (&)} )) A = ((Tr [F +i (0)F +j (& - £)] )) A (27) 

and the nuclear length is much larger than nucleon size due to confinement, -C we can approximate the 
quark-gluon correlation as [32l . [33j , 

„2 r 



4 

where 



2 2N C 

Afa N (*) 5j f I d&Mbf), (28) 



-Ja(pn, Zn)(N \...\N) = -Lpi(£ N )(N | • • • | JV>, (29) 



(27r)32p- J "^ IV,SIVM 1 17 2p 

denotes the ensemble (medium) average and pjy(£/v) is the spatial nucleon density inside the nucleus normalized to 
the atomic number A. The quark transport parameter <7f(£jv) is defined as 

9 2 JL,t^ [ d £~ ,an p ( (\\ i?g ( £—\ I art _ 2l ^ a ' n Art sr w i 



^ (6v) = ~2N C P ^ N) J t+ (N ' F +^ F +^~) I ^> = -^fpmN){xf 9 N (x)h=o, (30) 
and the gluon distribution function in a nucleon is 

<W = "/~ H+ e " P+r I ^(0)^ + CT (D I AO, (31) 

where the summation over the gluon's color index in the matrix element is implied in the definition of the gluon 
distribution function. 

The TMD quark distribution function in Eq. fT]) is a leading- twist result in terms of the momentum scale (Q 2 ) 
dependence of the DIS process. Among higher-twist corrections, one has neglected those from the transverse phase 
factors such as 



that the propagating quark accumulates in the above two-gluon correlation matrix element. They generally lead 
to contributions that are proportional to higher-twist nuclear parton matrix elements and are power suppressed 
(D(l/Q 2n ),n > 1. One can resurrect these higher- twist contributions by substituting k\ with its average value and 
setting the fractional momentum x — x± in the nucleon gluon distribution function in the quark transport parameter 
in Eq. ([50]) . This will introduce the Q 2 or energy dependence of the quark transport parameter (33[. For the rest of 
this paper we will focus our attention to the leading-twist TMD nuclear quark distribution. 

In the above approximation of the twist-four nuclear quark-gluon matrix we have neglected multiple-nucleon cor- 
relation in a large nucleus. Such an approximation is violated for small x where quark-gluon and gluon-gluon fusion 
from different nucleons become important and can lead to modification of the quark distribution function and gluon 
saturation in a large nucleus (34l . 13a . l36l | . We also neglected the real part of the nucleon gluon matrix elements which 
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is responsible for Pomeron-like elastic (or diffractive) scattering and the nuclear shadowing of the quark and gluon 
distribution functions [13, [M, H^, H(| One can effectivly take into account these effects by using a nuclear modified 
quark distribution function fq{xs) ^ Af^(xs) and saturated gluon distribution function in the transport parameter 
qp which could lead a non-trivial nuclear and energy dependence (33j . 

For other higher-twist nuclear matrix elements, we similarly consider only the case where quark and gluon fields 
are attached to different nucleons inside the nucleus. 



M 2 



Affi 



1 



(2n)W c 



((Tr W ± (y-)-V k± ))a 



2 1 1 



1 



(2ny.N, 



.((Tr 



Wi(0)-V fcj 



In 



(33) 



Again we will only keep the dominant terms that have the maximum nuclear enhancement for each given power 2n 
(or twist) of the transport operator. Such contributions come from contracting one pair of the gluonic fields with 
one nucleon inside the large nucleus. Because of color confinement, the relative longitudinal coordinate of the gluon 
pair is limited to the size of the nucleon while the average coordinate is set by the position of the nucleon which 
can be anywhere inside the nucleus. Therefore, each pair of the gluon fields will give rise to one power of nuclear 
enhancement factor Ra ~ A 1 / 3 . 

We will also neglect all terms that contain any power of the covariant derivative D±(0) in W±(0) since they are 
subleading in the nuclear enhancement comparing to the same twist nuclear matrix elements without any covariant 
derivatives. We call the above approximation maximal two-gluon correlation approximation since we reduce the 
multiple gluon correlations in the nucleus to products of two-gluon correlations that have the maximum nuclear size 
enhancement. The leading contribution to the 2n gluon correlation function is then 



— ((Tr 



In 



(2n)! 
2 n n\ 
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-1 



2N C 2p 



F + „(0)Fl(n I N) 



(34) 



This is essentially the extension of the approximation for twist-four quark-gluon correlation in a large nucleus to the 
case of quark-n-gluon correlation in which we assume the correlation of 2n gluon fields is approximately the product of 
n two-gluon correlators (or gluon distribution functions) . In the above equation, the combinatorial factor for grouping 
2n number of gluon field operators into n pairs, 



(2n-l)!! = 



(2n)! 



2 n n\ ' 

and the color factor for 2n gluon insertions, 



N C (N? - 1)» 



(iv c 2 - 1)" (2N c y ■ 



are used. Summation over polarization and color indices in the matrix elements for quark and gluon distributions are 
implied. For gluon propagation, the above color factor should be C\/(N^ — 1)™ instead. 

Using the definition of the quark transport parameter in nuclear matter <?_f(£at) as defined in Eq. (|30[) . we can now 
express the power expansion of the matrix elements as 



M 2n 



1 



(2n) 



dy - e ix P +y ( A i ^,(o)2_ 



de-f+x(r) ■ v» 



2n 



i>(y~) I A) 



(35) 



With the above simplification of the dominant contributions to the nuclear matrix elements, we obtain the transverse 
momentum distribution of the struck quark in DIS off a large nucleus from Eq. (|20| . 



tf(x,k ± )^ Aff(x) exp 



d£, N qFit 



S^ikA.) 



(36) 



in terms of collinear factorized (or transverse momentum integrated) quark distribution functions and the quark 
transport parameter qp which in turn is related to the gluon distribution density inside the nucleus. This result is 
also recently obtained by Majumder and Miiller [28[ via direct resummation of all twist diagrams in a covariant gauge 
calculation. 
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From Eq. (|36j) one can then calculate the total transverse momentum broadening of the struck quark due to multiple 
scattering inside the nuclear matter, 



1 



^2F 



Af q N (x) 



d 2 k ± klff(x,k ± ) = / dG&.(60> 



(37) 



which is the same as the twist-4 contribution [29[. Though other multiple parton scatterings contribute to the 
modified transverse momentum distribution, they do not affect the broadening of the mean transverse momentum 
squared within the maximal two-gluon correlation approximation. 
If we define transverse (coordinate) distribution as 



(38) 



The corresponding nuclear quark transverse distribution is, 

tfMx) = f ^e^ + ^(A|^(0)^e-^-^^)^(y-)M)«A/f(x)((e-^-^(^))) A 



exp 



J d£, N q F (£,N) 



(39) 



which has a Gaussian form in the transverse coordinate y± . It is then easy to obtain nuclear TMD quark distribution 
function as 



f q A (x, k x ) « Af?(x)-±— exp [-k 2 JA 2F ] , 
which is again a Gaussian with width given by the total transverse momentum broadening squared A2F- 



(40) 



B. Effect of nucleon TMD quark distribution 

In terms of twist expansion in the collinear factorization, one can consider the nuclear modified transverse momen- 
tum distribution in Eq. (|36p as the summation of all twist contributions. However, it contains only contributions 
with the dominant nuclear enhancement A"/ 3 in the 2(n+ l)-twist multiple-parton correlation inside a large nucleus. 
Such dominant multi-parton correlations in a large nucleus are shown to be made up of the products of leading twist 
nucleon parton distributions. We have neglected higher- twist contributions to the nucleon parton distribution, for 
example, the intrinsic transverse momentum of quarks inside a nucleon. 

Since the expression for the TMD nuclear parton distribution function in Eq. (|2T))) is general, it should contain these 
higher-twist effects. To isolate the contributions of the intrinsic quark transverse momentum inside a nucleon, we will 
make the following expansion of the matrix element in Eq. (f2"0")) , 



(A | $(p)l-e w ^y->^ VOT) \A)=J2 —i A I ^(°)^T \W±(y-) ■ V k± \ iP(y- 



H (2n)! 

n=0 y ' 



A) 



(2n)! 



^ (2n)! ^ (2m)!(2n - 2m)! N, 



1 r - -> i 2m 

((Tr W ± (y-)-V k± )) A 



xA(N | V>(0)-^- [Wl(y-)-V fcj 



1>(y-) | N) 



oo 1 

y — 



m=0 



_£1J_ 

2iV c 2^ 



V 



dZJrdtitifaXN | F +ff (0)F^(r) I N)^- 



E 



„+ r 



^(2n-2m)!^I^T 



^(y-) | AT) 



exp 



*;at<7f(£ 



A' 



A(jV | ^(0)^-e #J -(^> ffc -LV(y") I N), 



(41) 
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where the approximation for the expectation value of m pair of gluon fields in a nucleus in Eq. (|34p is used. Substitute 
the above matrix element into Eq. (|20p , we obtain the nuclear TMD parton distribution function, 



fq(x,k±) = Aexp 



(42) 



which is now related to the TMD parton distribution function of the nucleon, f^{x. k±). Comparing to the Eq. (|36p . 
in which the nucleon TMD quark distribution is assumed to be just a 5-function, the nucleon intrinsic transverse 
momentum in the above equation is the result of the inclusion of a subset of non-leading (in nuclear enhancement) 
higher-twist contributions. 

From the above nuclear TMD quark distribution one can derive a diffusion equation [HI, [2f| for the evolution of 
the quark transverse momentum distribution with the nuclear size (or propagation length), 



df q A (x,k A 



1 



9£ 



(43) 



JV 



with the diffusion constant given by the quark transport parameter <7f(£a0- Apparently, this is the reason why (If{£n) 
is often referred to as quark transport coefficient [21J . 

In coordinate space, nuclear quark transverse distribution can be obtained from Eqs. (|2"D|) and (f4"2"|) by partial 
integration, 



dy 



e lxp+y (A | V(0)— e" i#J - to 
2n 2 



i>{y-) | A) » Af^x^ae-^-^y-^ 



exp 



Af?(x,y ± ), 



(44) 



as the product of the nucleon transverse distribution and a Gaussian. The final quark transverse momentum distri- 
bution can then be obtained from the Fourier transform of the above, 



7tA 2 f 



d *e ±e -^-^ 2 ^f?(x,L), 



(45) 



as a convolution of the nucleon TMD quark distribution and a Gaussian with a width A 2 _f given by the path integral 
of the quark transport parameter or the total transverse momentum broadening squared [Eq. (|37|) ]. This is also a 
solution to the diffusion equation with an initial condition at £^ = given by the nucleon TMD quark distribution 
function, Af^(x,k±). 



C. Arbitrary gauge 

In an arbitrary gauge, one has to include both the longitudinal and transverse gauge links in the gauge invariant 
definition of the TMD parton distribution function in Eq. (fTJ). Following the same procedure as in the light-cone 
gauge, one can rewrite the nuclear TMD quark distribution function as 

ff(x,k ± ) = J ^|^d 2 ^e"^~-^e^ 

^^d 2 t ± e ix * + y-- ir ^(A | md ± )^e> S y^^£ TMD (0,y)^y-,y ± ) | A)S^(k ± ). (46) 

The transverse differentiation should act on both the quark field operator ip(y~ , y±) and the gauge link £tmd(0, y). 
We will use the following identity [Eq. (|118|) in Appendix A], 

id y± £ T MD(0,y) = £tmd(0, y)W±(y~ 7 y±) (47) 

with the gauge covariant form of the transport operator W±(y~, y±) given by, 

w ± ( y -,y ± ) = w ± (y-,y ± ) + g f dC£\(C,y-;y±.)F + ±(C,y±)£\\(C,y--,y±), (48) 
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which transforms like a covariant derivative D±(y~ ,y±) under any gauge transformation. One can recover from the 
above the equivalent identity in the case of light-cone gauge in Eq. (fig)) by setting A + = 0. These identities have 
been used to relate the T-odd and spin-dependent part of the quark distribution function to twist-3 parton matrix 
elements [U |42| . These twist-3 matrix elements are related to the first moments in transverse momentum of the 
TMD parton distribution function and the twist-3 contribution to DIS process. 

Completing the integration over the transverse momentum k± and coordinate y±, we can recast the nuclear TMD 
parton distribution function as, 



dy ixp+y' 

2ir 



(A | ^(0)^£||(0,ir)e^<*~>'^(ir) I A)5^{k 1 _). 



(49) 



It is in exactly the same form as the expression in the light-cone gauge in Eq. (j2"0|) except that there is the longitudinal 
gauge link which is necessary to ensure the gauge invariance of the above form of nuclear TMD parton distribution 
function. The transport operator Wj_(y~) = Wx(y~ , 0j_) is now given by its more general form in Eq. (|48[) . Integrating 
over the transverse momentum, one obtains the collinear factorized (or transverse momentum integrated) quark 
distribution function, 



dy 



(*) = / ^e^y-(A | m^w&y-My-) I a), 



(50) 



that is also gauge invariant under any arbitrary gauge transformation. In the explicit calculation of the transverse 
momentum broadening via cut-diagrams as in Ref. [28| , these gauge links arise from resummation of extra number of 
collinear gluons on either side of the cut in addition to the soft gluons with transverse momentum that contribute to 
the transverse momentum of the final quark. 

Following the same steps as in the case of light-cone gauge, we will be able to derive from Eq. (|49]) the nuclear 
modified transverse momentum distribution function in an arbitrary gauge as given in Eqs. (|42|1 and (|45|) . The 
corresponding quark transport parameter qf is simply replaced by a gauge invariant form 



2ir 2 a< 



2np~ i 



j xp +r {N i Tr[F +ff (0)£||(0,r)^(rK||(r,0)] | N), 



(51) 



F° a T a . Note that the above definition of 



where the gluon field is expressed in the fundamental representation F + , 
the gauge invariant gluon distribution function in the fundamental color representation is equivalent to the definition 
in the adjoint representation [43|, H|| , 



(N\F^mtaMr)F b + "(n\N) 



(52) 



where 



£|f(0,r) = exp 



-iff 



dCA c + (C,0 ± )t c A 



(53) 



with (t c A ) a i, = —ifabc, is the longitudinal gauge link in the adjoint representation. One can similarly introduce the 
transverse gauge link in the adjoint representation 



Cj_ (-oo;£_l,0 ± ) = exp 



-ig I d£±_ ■ A c ± (-oo,£ x )tA 



(54) 



and the TMD gluon distribution function, 

d 2 U d£ 



(2tt) 2 2ir P + e ' 1 ' 



(N | F^(0,0 ± )££ MDQb (0,£K ff (r,£±) | N) 



2np^ 



(N | F» (0) 



ab 



F b + °{C) | N)S^(k ± ) 



where 



C 



TMD 



Qb (o, o = rj^-oo, r ; Qx)^(-oo; a, ox)£f d6 (-oo, r ; a), 



'At 



(55) 



(56) 
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and 

wtbiCU) = iDt b (C,£±)+9 f dCcfl(c,r;^)Ft± ce (C,C±)cf eb (C,r;U), (57) 

J — oo 

is the transport operator in the adjoint representation. The covariant derivative and gluon field strength i^Tj_ in 
the adjoint representation are defined as 

Di ab (C,U) = S ab d ± + iglUC ,U)(t c A )ab, F^ ab = F° ± (t A ) ab , (58) 
respectively. The corresponding transverse coordinate gluon distribution is then 

Fl°(C)\N}. (59) 

ao 



2tTO+ X 



*+„(0) 



D. Effect of nucleon TMD gluon distribution 



In the maximal two-gluon correlation approximation, we approximate the higher-twist nuclear parton matrix el- 
ements with products of twist-two nucleon parton matrix elements [see Eqs. (j33| and (|34]) ] . Higher-twist nucleon 
quark matrix elements and quark-gluon correlations lead to intrinsic transverse momentum distribution inside the 
nucleon. However, we have so far neglected contributions from higher-twist nucleon gluon matrix elements that involve 
covariant derivative D± or multi-gluon correlation within a nucleon. These matrix elements have non-leading nuclear 
length dependence as compared to the products of twist-two gluon distributions. In order to consider the effect of 
these higher-twist nucleon gluon matrix elements, we separate the covariant derivative from the transport operator 
in light-cone gauge (to simplify notations), 



W ± (y-) = F ± (y-)+iD ± (y-); 
?±(V~) = 9 f dCF+±(C)- 

J —oo 

Note that [see Eq. (|120|) for arbitrary gauge in the Appendix] 



(60) 



D±i(v-)r±(y-) = g / drDUr)F+±(r)+^±(y-)D±.i(y-) = Di^ x (y-)+^(y-)D ±i (y- 



(61) 



where D^F = d±F + ig[A±,F] is the covariant derivative in the adjoint representation. Since one can factor out 
the covariant derivatives of the quark field together with other higher-twist nucleon quark-gluon matrix elements into 
the TMD nucleon quark distribution function, one can effectively replace the covariant derivative D± with its adjoint 
form in the the Taylor expansion of the nuclear gluon matrix element 



1 



f(v±) = Wc (P±e 



-iy ± -W±(y )\ 



(62) 



which should be a function of y\ because of the parity invariance of the unpolarized nucleus state. We will again 
neglect nucleon correlations and assume homogeneity in the nucleus. However, we now relax the maximal two-gluon 
correlation approximation to include higher-twist nucleon matrix elements that contain covariant derivatives and 
multiple gluon correlations. They are considered sub-leading in the nuclear length dependence in the above Taylor 
expansion. We denote fn(y±) as the nth term in the Taylor expansion of the matrix elements. Therefore, the linear 
term in y\ is 



y±-w ± (y~) )) A = 



vl 



1 



4 2N, 



■((WUy-)-wUy-)))A = --f^ F . 



(63) 



Note that the medium averaged value of matrix elements linear in J-±(y ) should vanish. 

For the quadratic term in y 2 ^, we first separate the gluon correlation into connected and disconnected parts, 



MvL) = 



i 



4\N r 



«Tr y x -W ± (y-) ))a=[^ 



1 



2\N r 



«Tr W ± {y~) ))a 



1 
2! 



A 2 



— ((Tr 



W x {y~) ))ac\=[^ 



—A 2 
2! 



2F 



A 



AF 



(64) 
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using the identity for generators of the fundamental representation, 

11 i 
T a T = -^-j-Sab + -^d abc T c + - fab c T c , 



(65) 



where the connected parts of the matrix elements ((■■ -))ac exclude the singlet contribution in the above color decom- 
position. We call twist-four quark transport parameter which contains all the twist-four nucleon gluon matrix 
elements in the connected part of the nuclear gluon matrix elements. 

In evaluating the connected parts of the nuclear gluon matrix elements, we will now adopt what we call extended 
two-gluon correlation approximation, in which we separate two gluon fields out of the nuclear matrix elements, 



MF 



^-((Tr [W ± (y-)] 4 )) AC » ^-((Tr [A(y~) ■ Wl(y-)A(y-)])) AC 
^((Tr {A • [t\ + 2f ± ■ iDi + (iDf ■ f ± ) + (*£>£ ) 2 ] Ta]))ac 
^T«n • [{F°)l b + 2(Wl ■ j%) ab + (iDf ■ j%) ab + {iDif ab ] F\_)) A 



2N, 



2 „A 



2np^ 



p A N ii 



N, 



2N r 



{N\F^)-[W°] 2 ab Fl^)\N). 



(66) 



Note that there are 2 pairs of gluons for the extended two-gluon correlation approximation and we have excluded the 
disconncted contribution (singlet) from the trace operation. We have defined the octect gluon field strength 



1 



•F. Lab — T^F\{dcab if cab) 

and the corresponding transport operator 

Again we have dropped terms linear in Tl{v~)- 

Similarly, one can also get the third term in the Taylor expansion of the transverse distribution, 



(67) 
(68) 



-l 



1 6 



((Tr ff ± -W ± (y~) ))a = - 



4 



3!A, 



1 ((Tr\wL(y-)] 6 )} A 



2 \ 3 

y± 



A 3 2F + 3A 2F — ((Tr [W±(y-)\ ))ac + jj-p* [W±(y-)\ ))ac 



4 



— A^ F + A 2 fA 4 j? + 7JjA 6F 



We again assume the extended two-gluon correlation approximation 

6 



(69) 



A 6F = 



^((Tr [W ± (y-)] 6 )) AC « ^((Tr [A(y~) ■ Wl(y-)fj_(y-) 



)ac= 2n} T1 



ah 



K))a 



(70) 



Following the same procedure and the extended two-gluon approximation, we can obtain other terms in the Taylor 
expansion and the final form of the transverse expansion, 



f(yi) 



i 



A 



2F 



el 

4 



2 r 



J vl 



vl 



— Alp + A 4F 



3 r 



— A 2F + A 2F A4 F 



1 

2!' 



^6F 



0(y\) 



exp 



(71) 



As one can observe, inclusion of higher-twist nucleon gluon matrix elements in the medium averaged products of 
the transport operator will give rise to a transverse-distance-dependent (TDD) quark transport parameter, 



A F (2/i) 



d^ N qF(^N,yi) 



d^dC. 



Trg 



2np+ 



2 n A tf \ n y 



(-i)" (vl 



2N n <—< n\ 

n=0 



.,2 \ n 



(N\FI ± (^).[W2] 2 2FI ± (^)\N) 



(72) 
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that depends on higher-twist nucleon gluon matrix elements. Such TDD quark transport parameter effectively con- 
tribute to power corrections to the nuclear transverse momentum distribution and render it a non-Gaussian form, 
especially in the small transverse momentum or large transverse coordinate region. However, these contributions 
have sub-leading nuclear length dependence as compared to those associated with the leading-twist quark transport 
parameter A2F- Since the leading twist (twist-2) transport parameter is proportional to the collinear nucleon gluon 
distribution function, such TTD transport parameter might be related to the TTD gluon distribution function. 
It is helpful, therefore, to make a similar Taylor expansion of the TDD gluon distribution function from Eq. (|59[) . 



2itp~ 



:(N\F% ± (0) 
' (-1)' 



F b + AO I AO 



E 



2np+ ^ (2n)! I ^x(o) • [y, • W?(n] ab FlAO I N) 
^tE (t)" {n 1 ^ (0) ' ^'tyiAD I N). 



(73) 



Comparing the above Taylor expansion of the TDD gluon distribution function to the TDD quark transport parameter 
in Eq. (f72")) . one can indeed relate the two, 



/4(&v)^[x/f(z,Z/i)]x=o, 



(74) 



if we approximate the octet gluon field strength with its adjoint value 



•Fxab — T^FA^abc ~ ifabc) ~ iFj_fa.bc — ^±ab- 



(75) 



We will refer to the above approximation and the extended two-gluon correlation approximation together as extended 
adjoint two-gluon correlation approximation. Under such approximation, we can resum contributions associated with 
the higher-twist nucleon gluon matrix elements to a TDD quark transport parameter which is related to TDD gluon 
distribution function. We would like to point out that under the extended adjoint two-gluon correlation approximation, 
one actully only includes a subset of higher-twist nucleon gluon matrix elements (see Appendix B). This is similar to 
the approximation we made in order to fact out the nucleon TMD quark distribution from the nuclear TMD quark 
distribution in Sec. Ill Bl 

Note that even though the higher-twist contributions we have considered so far lead to a non-Gaussian nuclear 
transverse momentum broadening, they do not contribute to the averaged transverse momentum broadening squared 
which is still given by the twist-2 transport parameter A2F- These higher- twist nucleon gluon matrx elements only 
contribute to higher moments of the transverse momentum broadening. 

One should also keep in mind that we have only considered the leading order in a s of the hard scattering. Higher- 
order contributions should also lead to leading-twist non-Gaussian components of the transverse momentum distribu- 
tion (|. 

Finally, in a finite nucleus, one has to take into account the finite number of nucleons A in the nucleus when 
factorizing the nuclear parton matrix elements into products of nucleon parton matrix elements. Such consideration 
will lead to a quark distribution function in the coordinate space (32j, 



Kvl) = i-(^-i) 



dt_N „/<t 2 \ 

—£-<i\£n,v±) 



~ te 2 



A-l 



(A-l)(A-2) 
2! 



—£-<U,N,y__) 



(76) 



which can be approximated as that in Eq. (|7ip for a large nucleus A > 1. 



III. QUARK PROPAGATION IN A THERMAL MEDIUM 



We can generalize our study of nuclear transverse momentum broadening to quark propagation in a hot medium 
such as the quark-gluon plasma produced in high-energy heavy-ion collisions. In this case, the initial quark production 
cross section is assumed to be factorized from the quark propagation in medium. 
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The thermal medium can be considered as an interacting gas of colored constituents with a local density Pn(£n), 
with TV now referring to the color constituents. The correlation length among these constituents are determined by 
the screening scale 1/p which could be longer than the inter-constituent distance as given by the temperature \/T in 
the weak coupling limit. Under such scenario, one can still apply the maximal two-gluon correlation approximation 
to the medium averaged multiple gluon matrix elements as we have used in the cold nuclear medium. 

With an initial condition of longitudinal momentum p + and zero transverse momentum, the final quark transverse 
momentum distribution can be easily read from Eq. (|42|) , 



= f(k x )5(x-l) 



exp 



5^(k ± ) 



7rA 



2F 



exp 



A 



2F 



or in terms of the transverse coordinate distribution 

,2 



/(£]_) = exp 



(77) 



(78) 



for the final quark. From the general form of the TMD quark distribution function in Eq. (|46j) , one can replace the 
nuclear state with a quark with momentum [p + , 0j_] and the medium A. Averaging over the initial state of both quark 
and the medium, one obtains the transverse coordinate distribution for a propagating quark as given by the medium 
expectation value of a pure gauge link, 



C\ (-oo, oo; 0±)£_l(-oo; 0j_, y±)C\\ (-oo, oo; y± 



(79) 



where we have assumed the quark is produced at oo and propagates toward — oo along the light-cone according to the 
convention used in this paper. 



A. Dipole model approximation 

In a covariant gauge (where the transverse gauge link becomes unity), the above becomes the Wilson line formulation 
of multiple scattering [23|, |4f| . Under a dipole model, the medium averaged Wilson line can be approximated [46[ as 



1 

At 



(Tr Ct (— oo, oo; 0±)£|| (— oo, oo; y±) 







) « exp 





(80) 



in terms of the dipole cross section <r(y±) and the medium density Pn(£,n)- Using the short distance form or the 
leading logarithmic approximation of the dipole cross section [20l | , 



Pjv(6v)o-(£i 



1 



one can obtain the expression in Eq. (|78|) for the transverse coordinate distribution f(y±). One can easily identify 
the first coefficient in the power expansion of the dipole expansion cif{£,n) with the quark transport parameter as we 
have defined in our twist expansion approach [Eq. (|51|) ]. Therefore, the maximal two-gluon correlation approximation 
for the dominant multiple gluon correlation in nuclear medium in our study here is equivalent to the dipole model 
approximation of the Wilson line approach [23|, |45[ when the short distance form of the dipole cross section is used. 

To relate our twist expansion result to the dipole model approximation beyond the maximal two-gluon correlation, 
we can use the identity in Eq. (| 1 18|) to recast the transverse coordinate distribution, 



e vx-d t±c t (_ QOj _ L )£ ± (-oo; ± , e±)£,| (-oo, oo; 



>^=o = i«Tre 



-iW±(oo)-y ± 



», (81) 



in terms of the transport operator Wj_(oo) [Eq. (|48p]. A Taylor expansion of the above distribution in y± in the 
extended adjoint two-gluon correlation approximation will lead to a transverse distribution as in Eq. (|7ip. 



f{y±)~ exp 



exp 



(-1)* 



(n-l) 



-A 



2nF 



(82) 
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Comparing the above distribution to the dipole model approximation in Eq. (|80[) . one can relate the dipole cross 
section to the nucleon TDD gluon distribution function, 



2 K 2 a s 



Va 



k/f(z,yi)]*=o. 



(83) 



This is exactly the cross section between a nucleon and a quark-anti-quark pair in a dipole configuration with transverse 
separation y± [H, E3] . 

Our calculation can also be extended to a gluon propagation. The results are the same and one only has to change 
the color factor to get the definition of the gluon transport parameter, 



QA(^N,y 2 ± ) 



N*-l 

Comparison to the approximation of the averaged Wilson line, 



1 



N? 



T^ Tr |£||(~ oo ' 0;0 -0 z: iK -oo ' 2/ ;£0j) ~ ex P 



(84) 



(85) 



in the study of gluon saturation in large nuclei [48| will relate the saturation scale Qsat(^i) with the path- integrated 
gluon transport parameter c}a(^n iV 2 ), 



QlAvi)= / d£ N q A ((N,y±) 



4:Tr 2 a s C A 
N? - 1 



(86) 



The transverse scale dependence of the transport parameter, the dipole cross section or the saturation scale in our 
calculation come from contributions of higher-twist nucleon gluon matrix elements and therefore are non-perturbative 
in nature. At very short transverse distance scale, radiative corrections will become important at the leading twist 
and they will give rise to a transverse scale dependence of the transport parameter that is governed by the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations 



ansport par 

SUM EH- 



B. Multiple gluon correlations in Af = 4 SYM 



In the Taylor expansion of the TDD quark transport parameter Ap(y 2 _^) in Eq. the coefficients are higher- twist 
nucleon gluon matrix elements which generally involve multi-gluon correlations of the medium. Therefore studying 
these power corrections to the transport parameter will shed light on multi-gluon correlations in a medium, especially 
a strongly coupled system when they become important. 

Non-perturbative calculation of the medium averaged Wilson line in Eq. (|81|) is difficult for a strongly coupled 
system. The developed technique of lattice QCD is not applicable because it is formulated in the Euclidean space 
and is only suited for the study of static thermodynamic observables. However, many transport coefficients, such 
as the shear viscosity to entropy density ratio r]/s [53[ and transverse momentum broadening of a heavy quark 
[54l [55| . have been studied for Af = 4 SYM theory in the large 't Hooft coupling (A = 3sym-^c) limit, employing the 
AdS/CFT correspondence [52 |. Though SYM is not exactly QCD, its study might provide indicative information on 
the properties of a strongly coupled system. 

Recently, the thermal averaged Wilson loop along the light-cone with longitudinal distance L~ and transverse 
separation y± was calculated [3fJ in the strong coupling limit of Af = 4 SYM theory. The corresponding transverse 
distribution from the calculated Wilson loop in the fundamental representation is 



f(y ± ) = cxp <^ -aVXL-T 



'irTyj 
2a 



1 - 1 



(87) 



where T is the temperature and a = y / 7rr(5/4)/r(3/4) w 1.311. Note a factor of l/\/2 is missing here because of our 
definition of light-cone variables. Expanding the exponent in a Taylor scries of the transverse distance y\, one has 



f(H±) = exp 



yl T _V\* 2 T* < 



2a 



E 

71=1 



Va 



n 2 T 2 
2a 2 



(-l)"- 1 |(2n-3)!!| 



(88) 
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Comparing to Eq. (|82)) . one can extract the leading- twist quark transport parameter in J\f = 4 SYM theory, 



-SYM 

If 



V\tt 2 T 3 
2a 



(89) 



which is half (versus 4/9 in QCD) of the gluon transport parameter as obtained in Ref. [301 ] . 

Comparing the power corrections, one can also extract higher-twist quark transport parameters in SYM (A/ - = 4), 



A 2rYf 



7T 2 T 2 

2a 2 



(2n-3)|! 



(n>l). 



(90) 



It is interesting to note that all higher-twist gluon matrix elements that include multi-gluon correlations are pro- 
portional to the leading-twist gluon matrix elements or two-gluon correlation. The coefficients are set only by the 
temperature of the medium and are independent of the coupling y / asYM an d number of colors -/V c . This implies that 
multiple gluon correlations in the strong coupling limit of SYM are as important as the two-gluon correlation. 

Phcnomcnologically, one can use the relationship between the transport parameter and the gluon distribution 
function in Eq. (jT4")l to obtain the TDD gluon distribution density in a Af — 4 SYM plasma, 



PN[xf g (x,yi)]™? 



8aN 2 T 



2a 



1 - 1 



(91) 



which is proportional to N 2 for fixed t' Hooft couling constant A. 



IV. SUMMARY AND DISCUSSION 



In this paper, we have derived a gauge invariant form of nuclear transverse momentum broadening distribution, 
utilizing the gauge invariant TMD quark distribution function from Ref. [lOj . We first express such TMD quark 
distribution function in terms of a sum of all higher twist and gauge invariant collinear parton matrix elements. 
These higher-twist parton matrix elements are expectation values of the moments of a transport operator W±(y) 
which generates the transverse momentum in a nucleon or nucleus when it acts on the parton field ip{y). The defined 
transport operator W±(y) transforms like a covariant derivative and the final expression is explicitly gauge invariant. 
With this general form of the TMD parton distribution function, one can then calculate any moments of the parton's 
transverse momentum in terms of the higher twist parton matrix elements. 

To calculate the nuclear broadening of transverse momentum distribution, we approximate the nuclear matrix ele- 
ments of n pair of parton field operators as a product of n nucleon parton distributions, neglecting nuclear bounding 
effect and multiple nucleon correlations which have non-leading nuclear size dependence in a large nuclear. In other 
words, multiple gluon correlations are assumed to be given as products of two-gluon correlations, which we called 
maximal two-gluon correlation approximation. With such approximated nuclear matrix elements that have the dom- 
inant nuclear size dependence of A n ^ 3 for fixed dimension of the multi-parton operators, we can express the final 
nuclear TMD quark distribution function in terms of the nucleon TMD quark distribution. This form also obeys 
a 2-D diffusion equation whose solution is a convolution of a Gaussian distribution function and the nucleon TMD 
quark distribution. The width of the Gaussian, or the mean total transverse momentum broadening squared, is just 
the path integral of the quark transport parameter qp, which is also defined in an explicitly gauge invariant form and 
is related to the local gluon distribution density. 

Under an extended adjoint two-gluon correlation approximation, one can resum some of the higher-twist nucleon 
gluon matrix elements to obtain a transverse-distance-dependent quark transverse parameter which is given by the 
TDD gluon distribution function. Such a TDD quark transport parameter will give rise to power corrections to the 
Gaussian form of nuclear transverse momentum distribution function. 

We compared our final results with that of the Wilson line approach to multiple parton scattering [23|, [4|| in the 
dipole model approximation. The two results are identical for short distance approximation of the dipole cross section, 
which is equivalent to the maximal two-gluon correlation approximation employed in our calculation. If we relax the 
maximal two-gluon correlation to include non-leading length-dependent contributions involving higher-twist gluon 
distribution functions, one can then relate the dipole cross section to the TDD gluon distribution function in the 
medium. 

We also compared our results with the AdS/CFT calculation [3fJ of the transverse distribution of a Wilson line for a 
propagating quark in the Af = 4 SYM theory, in particular the power corrections to the leading Gaussian distribution. 
We found that the SYM result indicates the importance of multiple gluon correlations in a strongly coupled system. 
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Though our final result for the nuclear modified transverse momentum distribution contains all higher-twist collinear 
nuclear parton matrix elements, it is still only the leading twist contribution in terms of power suppression 0(1/Q n ). 
One can follow the procedure as outlined in Ref. Q to compute higher-twist contributions to the momentum broad- 
ening which has the same nuclear length dependence but are power suppressed in the momentum scale of the hard 
processes as compared to the leading twist result obtained in this paper. 

The transverse momentum broadening we calculated in this paper is only valid in small transverse momentum 
region. At large transverse momentum radiative corrections become important and will lead to large logarithmic 
corrections to the transverse scale dependence of the transport parameter [33| . The form of nuclear broadening will 
also have large power corrections to the Gaussian form. 

The Gaussian form of the transverse momentum broadening in nuclear medium we discussed in this paper will 
also have phenomenological implications in the hard processes in hadron-nucleus and nucleus-nucleus collisions. Our 
study here justifies the phenomenological approach to the initial transverse momentum broadening [Blil . \El\ in p + A 
and A + A collisions in which a Gaussian form of the broadening is often used. Convoluted together with the power- 
law-like transverse momentum spectra due to hard scattering, the Gaussian broadening will naturally lead to Cronin 
enhancement of the final hadron spectra at small and intermediate transverse momentum. 
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APPENDIX A: Transverse derivative of a gauge link 

Consider a cap-like gauge link, 

C n = Clj_(-oc)C\\(-ao,y~; y± + dy±_)C d ±{y~) (92) 

with 

£d±{y~) = £±{y~;y± + dy±,yx,) = 1 - igdy± ■ Aj_(y~,y±) (93) 

and an infinitesimal transverse displacement dy± . One can break the longitudinal gauge link into a product of many 
small ones each with an infinitesimal length d£~ , 

£||(-oo,lT;£i_+d£i.) = •••£ d ||(i + l,i)---£ d ||(3,2)£ d ||(2,l), (94) 
£||(-oo,y-;yl) = • • • £„(* + 1,*) • • .£,,(3, 2)£,|(2, 1), (95) 

where Cm (i + 1, i) and C\\ (i + 1, i) are defined as 

£ d ,|(i + M) = C\\(& +1 ,^;y± + dy ± ), (96) 
£||(i + l,i) = £||(C+i>C;yU (97) 
and £7* = y~ — (i — l)d£~ . Inserting unit matrices 

i = £ dx (6)£||(i + M)£j(* + ( 98 ) 

between all neighboring links Cm + and Cm in £||(— oo, y~ ,y±+dy±) [Eq. as illustrated in Fig. 

except the last point where one instead inserts the unit matrix 

1 = £||(-co,-oo + dO£j | (-oo,-oo + ^ _ ) (99) 
after £^j_(— oo). One can then recast the cap-like gauge link, 

c n = £,| (-00,-oo + dr) •• + M)£ n (C)4i(M- i)£n(Ci) • • •£ n (£ 2 -)£||(2,i)£ n (£r). ( 10 °) 
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FIG. 2: Splitting the cap-like gauge link into chains of closed plaquette linked by short Wilson lines. 



as a product of closed plaquette, 

£□($-) = £j (t + 1, i)4 ± te+i)£ d || (t + 1, 

that are linked by short Wilson lines Cu (i + 1, i). 

Using the expansion of the closed plaquette [581 ] up to the linear term in d£,~dy±, 

one can expand £n up to the term linear in dy±, 

oo 
z=l 



= £||(-oo,j/ ;y^)-igdyj_ - / d£ £||(-oo,£ ,y_]_)£||(£ ,y ;y x )< 

J — oc 

Comparing to the direct expansion of £n in Eq. (|92|) in , 

£n = £||(-oo,y _ ; y ± ) +dyj_ ■ d y± C l \{-oo,y^;y±) 

- ig£\\(-°°,y~;y±)Aj_(y~ ,y±) + igA±(-oo,tf±)£\\(-oo,y~;jf±) , 

one obtains, 

d v ± £\\(~ oo, V~;yx) = -igA*j_(-oo,ifx)C\\(-ao,y~;y±) + C\\(-ao,y~;y±) 

igAj_(y~,y±) -ig [ dc ,y~;y±)F+±(C ,y±)£\\(C ,y~; y±) 



(101) 



(102) 



(103) 



(104) 



(105) 



One can also obtain the above transverse derivative via the direct expansion of the path-ordered longitudinal gauge 
link, 

d y± L\\ (-oo, y~;yx) = £H<7) n E / M >MG~, ia) ■ • ■ • • • yx) 

71=1 8=1 ^ 



EH5)"E / W^ttr.iTx)--- d + A ± (C,y±)-F + ^^y^) ...A+(c,yl), 

n=l i=l 



where 



50 r-oo 

[de]? = d^ d^ 

y Jy J y 



F +X (c.ta) = a + l ± (er,^)-a ± ^ + (er,?7±)- 



(106) 

(107) 
(108) 
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One can complete the integration of the terms with d + A±(£ i ,y±) by changing the order of integration, 



Jti+i 



-r 



for i = 1, 2, . . . , n — 1, with £ = — oo and 

rCi „ „ 

^„ d+A±(£_ n , y±) = A J _(^ n _ 1 ,y ± ) - A±(y , y±). 



(109) 



(110) 



We can rearrange the sum of terms associated with d + A±(£ >i ,y±) in the nth order (in the coupling g) of the expansion 



as, 



/ m^A + (^,y ± )...d + A ± (^,y ± )...A + (^,y ± ), 



n-l 

3=1 
n-l „_ 



n-2 



<<■ p — oo 
i=2 



^+(^„_i,j7±) [AL(£ n _i,£L)-AL(y 



i+l 



x A + (d, j?±) A ± (£"_ x , £0 - lj_ (C+i, £0 A + (C +1 , y ± ) / [^] ?+2 JJ A + (C7 , ft.) 



j=i+2 



+ / [«] 

= - / i^r 1 

1 /.-OO 

»=i j « 



A±(-oo,if±) - A ± (£ 2 ,y ± )] A+{Z 2 ,y ± ) J] A+(^ ,y ± ) 



J=3 



3 = 1 
i-1 



A±(y~,y±) + A±(-oo, $ ± ) / [dffi ]J A + (£r ,$ ± ) 

Jy- j=2 



3=1 



A + (£r,^),A ± (£r,^)l / [d^ J] ft.). (Ill) 

J y j=»+i 



The terms containing the commutator, , yj_), , can be combined with F + ±(£ i , yj_) in the (n — l)th 

order of the expansion to give the gluon field strength tensor, 



*+±(6 ,y±) = F +± (^ ,y 1 _)+ig[A + {^ i ,jf ± ),A±(£ i ,&.)]. 
Note that 

oo n—1 

EHs)"/ »+i II A + (^,y±) = C ll ^,y-,y±), 

71=1 J V j=i+l 



(112) 



(113) 



and 



i=i • / w 

OO /i QQ 

/ — oo 
d&C\{£i,-co,y ± ). 



^r^+(?r,y±)---^+fe-i,y±) 



(114) 
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One obtains now the transverse derivative of the longitudinal gauge link, 

d yx £\\(-ao,y-;ff±) = -ig f d^£\{C ,-co)F +± (C ,y±)£\\(^~ ,y~;jj±_) 

J —OQ 

+ig£\\(-oo,y~ ;y±)A±(y~ ,y±) - igA ± (-oo,y ± )£ n (-oo,y~;y ± ). (115) 
Another general form of the above identity is 

Dj_(yi,y±)£\\(yi,ir;y±_) = -ig f d4£ t l ^~,y 1 ~)F + _ L (^~,y±)£\\(Cy~;yi.) 

+£ ]l (y i T,y-;y ± )D x (y-,y ± ). (116) 
Using the above identity, the derivative operation on the gauge link in the TMD quark distribution function, 



(117) 



^tmd (0, y) = £^(-oo,0;6x)£±(-oo;y±,^±)£^(-oo,y ;y±), 
will yield, 

d y± £ T M-D(0,y) = £tmv(0, y)d y± + £^(-<x,0;0±)£]_(-oo;yx,0x)d y± £\\(-cx} ) y~ ;y±) 

+ig£\(-oo, 0; ± )£ 1 L (-cx); y±, x )A±(-oo, y±)£\\ (-oo, y~ ; y±) 

= £ T md(0,?/) D±(y~,y±)-ig d£~ £\(C ,y~ ;y±)F +± (C ,y±)£\\(C ,y~;y±_) ■ 

J — oo 

If we define 

F±(y ,y*±) = g d£ £](( ,v ;y±)F+±(£ ,y±)£\\(€ ,y ;v*±), 

J —oo 

it is easy to use Eq. (|116[) to show 

Dxi(y,y±)A(y-,yx)=f±(y-,yx)D X i(y,y±)+9 f dr£\(r,y'-,y±)DiiF +± (r,y±)£\\(r,y--,yx), (120) 

where D±F = d±_ + ig[A±, F) is the covariant derivative in the adjoint representation. 



(118) 



(119) 



APPDENIX B: Extended two-gluon correlation approximation 

In this Appendix we will examinchighcr-twist nucleoli gluon matrix elements that are neglected in the extended 
two-gluon correlation approximation in Sec. Ill Dl 

In the connected part of the twist-four gluon matrix elements one can extend the product of transport operator 
without assuming extended two-gluon correlation approximation, 

A 4 f = ^((Ti-[W ± (y-)] 4 ))Ac 



(121) 



Terms linear in J r ±(y ) are dropped since they vanish after the medium average. The following identities along 
medium averaged matrix elements are also used, 



«Tr T x ■ {DlYT x )) A = <(Tr (Df)^ ± ■ F ± )) A = -«Tr {D± ■ T x ) ■ [D± ■ T ± ) )) 



(122) 
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«Tr 



((Tr Tx. ■ (iDi ■ T ± )f ± )} A = ((Tr 



(123) 



One can see that there are differences between the higher-twist gluon matrix elements in the above expansion and that 
under extended two-gluon correlation approximation in Sec lIIDl which include some extra higher-twist gluon matrix 
elements. 
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